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Abstract This paper presents a structural model called 
“PB3D” to perform nonlinear time history analyses of post 
and beam timber buildings under seismic loads. The model 
treats the three-dimensional structure as an assembly of 
roof/floor diaphragms and wall subsystems. The roof/floor 
diaphragms are modeled by beam elements and diagonal 
brace elements in order to take into account the in-plane 
stiffness. The wall system consists of vertical beam ele¬ 
ments, for wall posts, as well as nonlinear shear springs to 
consider the contribution of diagonal wall bracing mem¬ 
bers or sheathing panels. The hysteretic characteristics of 
the shear springs are represented by a simplified, 
mechanics-based model named a “pseudo-nail”. Standard 
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finite element procedure is used to construct the system’s 
equation of motion, which is solved by Newmark’s inte¬ 
gration. The model was verified against shake test results of 
a three-story post and beam building subjected to strong 
ground motions scaled to the 1995 Kobe earthquake. 
Model predictions agreed very well with the test results in 
terms of base shear forces and inter-story drift responses. 
This model provides a robust and efficient tool to study the 
seismic behavior of post and beam timber buildings. 

Keywords Seismic performance • Timber buildings • 
Nonlinear analysis • Post and beam construction 

Introduction 

The 1995 Hyogo-ken Nanbu (Kobe) earthquake raised a lot 
of concerns about the seismic safety of traditional post and 
beam (P&B) timber construction in Japan. The poor per¬ 
formance of the severely damaged P&B houses can be 
attributed to some structural design issues such as heavy tile 
roofs, weak first stories, irregular shear wall layout, and 
inadequate provision of inter-story and foundation anchor¬ 
age. Poor maintenance was also a culprit [1]. Since the Kobe 
earthquake, a lot of efforts have been directed to study the 
seismic behavior of the P&B buildings and the Building 
Standard Law (BSL) in Japan has also been upgraded. 
Similar to the wood-frame construction in North America, 
P&B timber buildings in Japan are also box-type structures 
consisting of two-dimensional horizontal and vertical 
assemblies such as walls, floors, ceilings and roofs. Nail 
fasteners are also extensively used to attach sheathing panels 
to timber frames. However, the timber frames of the P&B 
buildings are usually constructed by members with rela¬ 
tively large cross sections, for example, 105 mm x 105 mm 
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for wall posts and 105 mm x 210 mm for floor/roof beams. 
A special feature of the P&B construction is the extensive 
use of traditional mortise-and-tenon joinery reinforced by 
metal hardware. Thus, it is very challenging to develop 
detailed finite element (FE) models to simulate the seismic 
response of such complicated building systems. 

Performance-based seismic design of wood buildings 
requires reliable and efficient numerical models to perform 
nonlinear time history analyses. Compared with the state- 
of-art numerical models for steel or concrete structures, 
very limited numerical models have been developed to 
study wood structures. A simple model to study the seismic 
responses of a wood building might be a multidegree-of- 
freedom (MDOF) model in which each floor or roof is 
lumped into a structural mass point and nonlinear shear 
springs are used to connect floors and roof. Only horizontal 
degrees of freedom (DOF) are considered. The character¬ 
istics of each shear spring are calibrated by the total con¬ 
tribution of the shear walls at each story. Agawa and 
Miyazawa [2] used such a MDOF model to study existing 
residential buildings in Japan. Although MDOF models are 
very computationally efficient for nonlinear time history 
analyses, they are highly simplified and cannot consider 
torsional effects and vertical effects in the buildings. 
Researchers also developed models which are able to 
consider the hysteresis of individual shear walls. Folz and 
Filiatrault [3] reported a benchmark modeling study of a 
two-story building tested in the CUREE-Caltech wood- 
frame project. Blind predictions with the assumption of 
rigid roof/roof diaphragms were reported. One of models 
was built in the RUAUMOKO program [4] and two were 
built in DRAIN-2DX and DRAIN-3D program [5]. 
Recently, van de Lindt et al. [6] reported a model called 
SAPWood to perform time history analyses of wood-frame 
buildings. In this model, the assumption of rigid floor/roof 
diaphragms was also used. The model consists of nonlinear 
shear springs for shear walls and non-symmetric linear 
springs to provide vertical restraints for the buildings. The 
assumption of rigid roof/floor diaphragms can greatly 
simplify the analyses since only three DOFs (two transla¬ 
tional and one rotational) need to be considered for each 
floor or roof. The model accuracy greatly relies on the level 
of the accuracy of shear wall hystereses implemented into 
the building model. 

Detailed FE models have also been developed to model 
wood buildings. He et al. [7] developed a sophisticated 
model called LightFrame3D for wood-frame buildings. In 
this model, a wood-frame building was modeled as an 
assembly of generic super-elements consisting of framing 
elements, panel elements and nail connection lines. Nak- 
agawa and Ohta [8, 9] and Nakagawa et al. [10] developed 
a model to simulate the collapsing process of traditional 
Japanese wood houses under severe ground shaking using 


the extended distinct element method since the authors 
believe that it is difficult to simulate such an ultimate 
failure process using common finite element methods. This 
model provides a very insightful tool to study the rela¬ 
tionship between seismic intensity and collapse limit of 
timber buildings. Commercial structural analysis software 
has also been used to model wood buildings, although most 
of them might have limited ability to represent the hys- 
teretic characteristics of wood connections/walls, such as 
strength/stiffness degradation and pinching effect. Masa- 
lam et al. [11] developed a dynamic model in SAP2000 for 
a three-story wood-frame building tested in the CUREE- 
Caltech wood-frame project. Collins et al. [12] developed a 
static model in ANSYS for a one-story wood-frame 
building. Lam et al. [13] also developed a static model in 
ANSYS for a one-story P&B building. Xu and Dolan [14] 
developed a dynamic FE model in ABAQUS to study the 
seismic behavior of a two-story wood-frame building in 
which the modified Bouc-Wen-Barber-Wen (BWBN) 
model was integrated to represent the hystereses of shear 
walls. 

Typically, detailed FE models are very computational 
intensive to perform nonlinear time history analyses of 
complicated wood systems although they are more com¬ 
prehensive and more accurate to predict structural 
responses. Therefore, they may not be suitable for seismic 
reliability evaluations which require a robust and efficient 
computer model to run a large number of simulations, 
considering the uncertainties in ground motions and 
structural capacity. In this study, a computer model called 
“PB3D” is developed to model the seismic behavior of 
the P&B timber buildings. This model aims to capture the 
global characteristics in the seismic response of the 
buildings with reasonable accuracy and computational 
efficiency by incorporating simplified structural component 
models for the roof/floor diaphragms and shear walls. For 
model verification purpose, shake table test results of a 
full-scale three-story building, excited by strong ground 
motions, were used. 

Theory 

Figure 1 illustrates the “PB3D” model schematics of a 
simple two-story P&B building. The floor and roof dia¬ 
phragms are modeled by 3D linear elastic beam elements 
and diagonal braces which are calibrated truss elements to 
maintain the same in-plane stiffness of diaphragms in real 
buildings. Shear walls are modeled by vertical beam ele¬ 
ments and nonlinear shear springs. The nonlinear springs 
are to represent the lateral resistance provided by wall 
bracing elements such as diagonal bracing members, 
structural sheathings (i.e., plywood and oriented strand 
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boards/OSB), and gypsum wallboards (GWB). The load- 
deformation hysteresis F(A) of the shear springs are 
calculated by a mechanics-based subroutine named “pseudo- 
nail”. 

Shear walls 

In many building models, shear walls are represented by 
shear springs with F( A) being the load-drift relationship 
and defined by mathematical functions. These functions 
associated with loading and unloading rules are usually 
obtained by empirically fitting hysteretic loops of shear 
walls from cyclic testing or detailed shear wall modeling. 
Numerous empirical hysteresis models for wood shear 
walls have been discussed in literature [15-18]. In the 
“PB3D” model, the shear walls are modeled by two types 
of elements: beam elements for the vertical posts to carry 
vertical loads and provide vertical restraints; and nonlinear 


shear springs for lateral resistance. The P — A effect is not 
considered. Since the contribution from wall posts to lateral 
resistance is small, the lateral resistance of the building is 
primarily provided by the shear springs. The uplifting 
characteristics of the building are not modeled by adding 
linear/nonlinear vertical spring elements. Instead, the 
uplifting is simply prevented by wall post elements which 
are fully end restrained onto the foundation or stories. The 
hysteresis F(A) for each shear spring, i.e., each shear wall, 
is calculated by the “pseudo-nail” model, as shown in 
Fig. 2. The concept of this model was motivated by the fact 
that a nailed connection and a wood shear wall demonstrate 
similar hysteretic features such as strength/stiffness deg¬ 
radation and pinching effect. The algorithm for the 
“pseudo-nail” model is a revised version of a nail con¬ 
nection model called HYST [19], originally developed to 
model wood connections with metal fasteners, except that 
the “pseudo-nail” model parameters need to be calibrated 
to represent a shear wall hysteresis instead of a nail hys¬ 
teresis. These parameters include the nail length L and 
diameter D and six parameters to describe the compressive 
properties the surrounding embedment medium. The 
“pseudo-nail” model is very computationally efficient and 
is believed to be more robust than the empirical curve-fitted 
models because, to some extent, it embodies the internal 
mechanism of a wood shear wall to resist lateral loads such 
as the interaction between wood medium and metal fas¬ 
teners. Furthermore, this model can adapt to any loading 
protocols while using a fitted hysteresis for any loading 
implies an extrapolation from the protocol used for the 
fitting. The details of the “pseudo-nail” model have been 
discussed elsewhere [20, 21]. Recently, Li et al. [22] 
modified the original HYST algorithm by adding the 
stiffness degradation in wood embedment medium under 
cyclic loading in order to fully address the strength/stiff¬ 
ness degradation and pinching effect in wood connections. 
The modified algorithm has also been applied into the 
“pseudo-nail” algorithm implemented into the “PB3D” 
model. 


Fig. 2 A shear wall hysteresis 
represented by “pseudo-nail” 
model 
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Roof/floor diaphragm 

Observations from past earthquake experiences and shake 
table tests show that when wood buildings are excited by 
strong ground motions, roof and floor diaphragms usually 
do not experience high nonlinearity. Severe structural 
damage and most of the input energy dissipation take place 
in the wall systems. Therefore, in the “PB3D” model, roof/ 
floor components are modeled by linear elastic 3D beam 
and truss elements. The roof and floor beams with large 
cross sections are modeled by the beam elements with each 
node having three translational DOFs and three rotational 
DOFs. Other components such as joists/rafters, sheathing 
panels and nail connections are condensed and simply 
modeled by diagonal truss elements with each node having 
three translational DOFs. The connections between the 
beam elements and the truss elements are simply assumed 
to be rigid. The cross sections of truss elements are cali¬ 
brated to match the in-plane stiffness of the floor/roof 
diaphragms in actual buildings. The element stiffness 
matrices can then be calculated and assembled into a global 
system stiffness matrix. 

Formulation of system equation 

The system equation of motion for the “PB3D” model is 
formulated as: 

Mil + Cu +/ s (m, u) = — Mu g (1) 

where M and C are, respectively, the mass and the damping 
matrix; u , u and il are the displacement velocity and 
acceleration vector of the finite element assemblage; il g is 
the ground acceleration; / s (w, u ) is the resistant force vector 
depending on the load history of the nonlinear system. 
Damping mechanism in wood buildings is very 
complicated. In general, under strong ground motions, 
the major energy dissipation is attributed to structural 
nonlinear deformations while the secondary energy 
dissipation mechanism can be approximately accounted 
by the viscous damping. In this model, mass-proportional 
and stiffness-proportional damping are used to account for 
the viscous damping: 

C = otM + /]K (2) 

where M is the lumped mass matrix and K is the stiffness 
matrix. The mass-proportional damping ratio a and stiff¬ 
ness-proportional damping ratio / can be determined based 
on test results. 

Newmark’s integration with constant acceleration 
assumption for each time step is used to calculate the time 
history response of the system. Given a short time interval 
At, the incremental form of Eq. (1) is 


MAilt + CAili + (A/^. = A Pi (3) 

(A/); can be approximated by tangent stiffness matrix 
of the system as 

(A/s),. = (Ki) T Aui (4) 

in which (/L)t is the global tangent stiffness matrix which 
is an assembly of element tangent stiffness matrices fol¬ 
lowing standard finite element procedures. 

In the “PB3D” model, for computation efficiency, the 
stiffness matrices for the linear elastic elements are 
established in a pre-process procedure. Thus, during iter¬ 
ations, only the stiffness matrices of the nonlinear elements 
need to be updated for each time step. Newton-Raphson 
iterations are used to solve for the incremental nodal dis¬ 
placement vector. Both energy and force convergence 
criteria are used to assure a converged solution. 

Model verification 

Shake table test results of a three-story building were used 
to verify the “PB3D” model. The building was tested on an 
8 m x 8 m large-scale shake table in the Public Works 
Research Institute (PWRI) in Japan. The building was 
designed according to allowable stress design method for 
wood post and beam construction [23] and the total 
effective lengths of shear walls as well as the eccentricity 
ratios also satisfied the Building Standard Law (BSL) in 
Japan. Figure 3 shows the building on the shake table and 
the “PB3D” model snapshot. This building had a plan size 
of 7.28 m x 7.28 m and was 9.5 m high. The timber frame 
was constructed by Canadian Hem-fir timbers with grade 
stamped as Canada Tsuga El20. The shear wall system, as 
shown in Fig. 4, consisted of OSB-sheathed walls, cross- 
braced walls, GWB-sheathed walls as well as their com¬ 
binations to enhance the lateral capacity. JAS Grade 4 OSB 
panels in thickness of 9 mm were used as wall sheathings 
and fastened with the timber frame with JIS CN50 nails 
(50.8 mm in length and 2.87 mm in diameter) spaced at 
100 mm. JIS GB-R boards in thickness of 12.5 mm were 
used to provide fire resistance as well as additional lateral 
resistance. The cross braces were 90 mm x 45 mm timber 
members with ends connected by EG gussets (equivalent to 
BP-2 plates). S-HD25 hold-downs were used to restrain the 
walls from uplifting. According to the wind loads, the 
required minimum total effective length of shear walls in 
the first story by the BSL was 31.0 m along the v or 
y direction. However, according to the seismic loads, the 
required length of shear walls in the first story was only 
17.9 m along v or y direction. Therefore, the design of 
shear walls in this building was actually governed by the 
wind loads. In this building, the total effective length of 
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Fig. 3 Three-story building and the “PB3D” model 



1 st story 


— OSB-sheathed wall; 

Sujikai wall (diagonal braced); 

— GWB-sheathed wall. 




3 rt * story 


Fig. 4 Shear wall layout of the building 


existing shear walls in the first story, excluding the walls 
with openings, was about 71.9 m along the x direction and 
78.4 m along the y direction, ratios of 2.32 and 2.53 of the 
code minimum by the wind loads, and ratios of 4.02 and 
4.38 of the minimum requirement stipulated by the BSL. 

This building also had a wide carport opening in the first 
story which greatly increased the building eccentricity and 
caused significant torsional effect under seismic loading. 

In the “PB3D” model, the diagonal truss elements for 
the floor diaphragm were calibrated by a static test of a 
floor segment with a plan size of 3.6 m x 5.4 m (Fig. 5). 

All the floor beams and joists were Canadian Hem-fir 
members. Plywood panels with thickness of 16 mm were Fig. 5 In-plane static test of a floor diaphragm segment 
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sheathed onto the frame by JIS CN50 nails spaced at 
150 mm on edge and 300 mm on field. Third point loading 
was applied via a steel loading beam onto one floor side. 
Both supporting sides were tied down on the strong con¬ 
crete floor in the lab. The in-plane stiffness of the floor 
segment along both horizontal v and y directions was tes¬ 
ted. However, along the deeper v direction, considering the 
floor aspect ratio, the concept of treating this floor segment 
as a “deep” beam may not be appropriate since the loading 
effect was highly concentrated in an area close to the 
loading points. The test results indicated that the in-plane 
stiffness of the floor segment was about 8.45 kN/mm along 
the y direction. Accordingly, the diagonal bracing elements 
of floor diaphragms in the building model were calibrated. 

The in-plane stiffness of a typical roof segment in the 
P&B construction was estimated by detailed FE modeling 
using the commercial software ANSYS, as shown in Fig. 6. 
The modeled roof was also constructed by Canadian Hem- 
fir members and sheathed by 12.5 mm thick plywood 


panels. JIS N50 nails were spaced at 100 mm on edge and 
200 mm on field. Beam elements and panel elements 
and link elements were used to model the roof beams 
and rafters, sheathing panels and nailed connections. The 
cross section of the diagonal truss elements was calibrated 
by matching the in-plane stiffness of these two structural 
systems. 

The “pseudo-nail” model parameters for individual 
shear walls were calibrated by a test database of P&B shear 
walls constructed with different wood species, sheathing 
materials, nailing schedules, and construction methods. In 
such tests, the load-drift backbone curves from reversed 
cyclic tests are usually used to determine the shear wall 
design values considering yielding, ductility and energy 
dissipation and deformation limit of the walls. The back¬ 
bone curves of the load-drift loops were used to calibrate 
the “pseudo-nail” model parameters which were imple¬ 
mented into the building model. Table 1 shows the cali¬ 
brated “pseudo-nail” model parameters for the six types of 
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Fig. 6 Roof diagonal bracing calibration by FE models 


Table 1 “Pseudo-nail” model parameters in three-story building 


Model Parameters 

0.91 m 

X-brace + 

OSB + GWB 

0.91 m 

X-brace + 
GWB 

0.91 m 

Doubl. 

GWB 

1.82 m 

Doubl. 

GWB 

0.91 m 

OSB + 

GWB 

1.82 m 
OSB + 
GWB 

Embedment 

Qo (kN/mm) 

2.38 

1.06 

0.25 

0.5 

1.2 

2.4 

Qi (kN/mm 2 ) 

0.025 

0.02 

0.002 

0.006 

0.006 

0.01 

& 

1.18 

1.18 

1.6 

1.6 

1.3 

1.3 

K (kN/mm 2 ) 

0.28 

0.2 

0.05 

0.1 

0.2 

0.3 

£> max (mm) 

75 

65 

70 

70 

75 

75 

Nail (mm) 

L 

300 

275 

225 

250 

275 

300 

D 

6.0 

5.5 

4.5 

5.0 

5.5 

6.0 
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full-height walls used in this building. The definitions of 
these parameters have been given in literature [20, 21]. 
Figure 7 shows the model predicted load-drift hystereses 
subject to a standard cyclic loading protocol. The hysteretic 
characteristics of strength/stiffness degradation and pinch¬ 
ing effect have been full addressed. 

In this study, the building model did not consider the 
contributions from exterior wall claddings (FMQ5MIV 
siding panels manufactured by AT-WALL) to the building 
lateral resistance because these claddings were simply 
hanging on the exterior walls by metal clips and it is 
believed that their contributions to the shear wall resistance 
are very small compared with those from the OSB shea¬ 
things, cross braces and GWBs. It should be noted that the 
building model also did not consider the contributions from 
the walls with openings. One reason is that this building 


had a limited number of walls with window/door openings. 
For example, using the opening reduction factor K 0 stipu¬ 
lated by the guideline from Seismic Evaluation and Retrofit 
for Wood Houses by Japan Building Disaster Prevention 
Association [24], the calculated effective length of walls 
with openings in the first story was 3.185 m along the 
v direction and 4.55 m along the y direction, approximately 
4.4 and 5.8% of the total effective length of full-height 
walls with OSB sheathing, cross braces, and GWBs along 
the corresponding direction in the first story. Therefore, the 
lateral resistance of this building was mainly governed by 
the full-height shear walls. Of course, the model predic¬ 
tions of the seismic response of this building would be 
more accurate should the actual shear resistances of these 
walls with openings be available and considered by the 
model as well. In other situations, if a building has 
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Fig. 7 “ Pseudo-nail” model predicted load-drift hystereses 
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Fig. 8 Major failure modes in 
shear walls 



150% Kobe JMA, X direction 


200% Kobe JMA, X direction 




150% Kobe IMA, Y direction 


200% Kobe JMA, Y direction 




Fig. 9 Base shear versus 1st story drift responses under 150 and 200% Kobe JMA shaking 
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significant amount of walls with openings or the wall fin¬ 
ishing materials have significant contributions to shear wall 
resistance, the seismic simulations should take into account 
their influence on the seismic response of the entire 
building. 

The total weight of the building was about 341 kN. The 
structural mass (including the structural self-weight and 
additional mass to simulate live loads) was 228 kg/nT for 
the 2nd floor, 225 kg/m 2 for the 3rd floor and 159 kg/m 2 
for the roof. The input ground motions were the three axial 
1995 Kobe JMA ground motion records including EW, NS 
and UD components. The weak v direction (due to the 
carport opening) of the building was excited by the NS 
component which has a higher PGA than the EW compo¬ 
nent. The building was tested with 150 and 200% Kobe 


JMA records, respectively. The recorded PGAs were 1.35 
and 1.75 g along the v direction of the building. During the 
150% Kobe JMA shaking, the peak inter-story drift (PID) 
of the entire building occurred at the carport opening in the 
first story and was about 59 mm, a drift ratio of 2.1%. No 
serious structural damage was observed except some corner 
breakage of GWBs, minor damage at OSB comers, and 
some detachment of the stucco siding in the first story. No 
repair was thus made after this shaking. However, during 
the 200% Kobe JMA shaking, at the carport opening, the 
PID reached 196 mm, a very high drift ratio of 7.2%. This 
building was at a near collapse state. Most of the shear 
walls in the first story experienced severe structural dam¬ 
age. As shown in Fig. 8, the OSB sheathing was detached 
from the frame due to the complete failure of nail lines and 
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the diagonal bracing members were broken due to the out- 
of-plane buckling failure. 

In the “PB3D” model simulation, structural mass was 
lumped onto roof/floor nodes based on their tributary areas. 
Combined ground motions of the 150 and 200% Kobe JMA 
shakings were input as one time history analysis since no 
structural repair was made between these two shake table 
tests. A damping ratio of 5% was assumed in the model 
simulation based on a previous free vibration test on a one- 
story post and beam building in which a damping ratio of 
4.7% was recorded [20]. The 5% damped model predic¬ 
tions also indicated that the PID of the entire building 
occurred at the carport opening side in the first story. 
Figure 9 shows the comparison between the model pre¬ 
dictions and the test results in terms of the relationship 
between the base shear and the first story drift. In general, 
good agreements can be observed. During the 150% Kobe 
JMA shaking, the model predicted maximum base shear, 
along the v direction, was 434 kN and the prediction error 
was only 0.5%. During the 200% Kobe JMA shaking, the 
model predicted maximum base shear, along the v direc¬ 
tion, was 544 kN and the prediction error was about 
—3.5%. The model predicted inter-story drift response also 
agreed very well with the test results, especially for the 
150% Kobe JMA shaking. Figure 10 shows the predicted 
first story responses at different locations where displace¬ 
ment sensors were installed compared with the test results 
from the 150% Kobe JMA shaking. The predicted PID at 
the carport opening side was 56.6 mm, 4.0% lower than the 
test results. Figure 11 shows the predicted first story drift 
responses along the v direction compared with the test 
results from the 200% Kobe JMA shaking. The predicted 
PID was 194.5 mm, 0.7% lower than the test results. 
Table 2 gives the model predicted PIDs at different 
building locations compared with the tests results. 





Fig. 11 Inter-story drift responses under 200% Kobe JMA shaking 


Table 2 Model predicted peak 
inter-story drifts versus test 
results 


Shake table test 

Sensor 

location 

Peak inter-story drift (mm) 



1st story 

2nd story 

3rd story 

150% Kobe JMA 

XI 

56.6 (-4.0%) 

25.6 (18.0%) 

11.3 (-28.5%) 


X2 

34.0 (-12.6) 

23.4 (0.4%) 

9.9 (-2.0%) 


X3 

24.9 (15.7%) 

18.7 (12.5%) 

7.9 (36.3%) 


Y1 

28.7 (28.4%) 

11.9 (-1.6%) 

3.9 (17.6%) 


Y2 

20.5 (24.3%) 

18.9 (19.6%) 

5.2 (27%) 


Y3 

22.4 (11.1%) 

23.0 (24.6%) 

6.3 (-11.8%) 

200% Kobe JMA 

XI 

194.5 (-0.7%) 

87.5 (-4.4%) 

24.0 (-34.4%) 


X2 

131.0 (1.9%) 

64.4 (-8.7%) 

21.9 (-4.8%) 


X3 

74.3 (-7.7%) 

29.7 (-14.3%) 

17.5 (23.2%) 


Y1 

51.8 (17.2%) 

34.8 (7%) 

7.0 (22.8%) 


Y2 

50.8 (23.0%) 

31.2 (18.0%) 

9.5 (3.2%) 


Y3 

48.4 (13.9%) 

25.9 (9.7%) 

7.3 (12.3%) 
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The numbers in parentheses give the prediction errors in 
percentage. In general, model predictions agreed very well 
with the test results. 

Conclusions 

This paper presents a computer model called “PB3D” to 
simulate the seismic response of post and beam wood 
buildings commonly used in Japan. The roof/floor dia¬ 
phragms are modeled by beam elements and truss elements 
considering the in-plane stiffness of diaphragms. The shear 
walls are modeled by vertical beam elements as wall posts 
and nonlinear shear springs. The hysteresis of the nonlinear 
shear springs is represented by a mechanics-based algo¬ 
rithm called a “pseudo-nail” model which can be cali¬ 
brated by shear wall test results or detailed wall models. 
Model predictions of a full-scale three-story building were 
compared against the shake table test results. The predic¬ 
tions agreed very well with the test results, especially in 
terms of the base shear forces and the first story inter-story 
drift responses. The torsional effect was also captured very 
well by the model. 

For seismic reliability analysis of timber buildings, a 
major computational demand is to establish a seismic 
response database, considering different seismic events and 
structural characteristics. The “PB3D” model presented in 
this paper provides a robust and efficient tool to estimate 
the seismic response of the post and beam timber buildings, 
facilitating the development of reliability-based assessment 
of the seismic safety of such building systems. 
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